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We study spin relaxation and dynamics of collective spin excitations in correlated double-exchange 
ferromagnets. For this, we introduce an expansion of the Green's functions equations of motion that 
treats non-perturbativerly all correlations between a given number of spin and charge excitations 
and becomes exact within a sub-space of states. Our method treats relaxation beyond Fermi's 
Golden Rule while recovering previous variational results for the spin-wave dispersion. We find 
that the momentum dependence of the spin-wave dephasing rate changes qualitatively due to the 
on-site Coulomb interaction, in a way that resembles experiment, and depends on its interplay 
with the magnetic exchange interaction and itinerant spin lifetime. We show that the collective 
spin relaxation and its dependence on the carrier concentration depends sensitively on three-body 
correlations between a spin excitation and a Fermi sea electron and hole. The above spin dynamics 
can be controlled via the itinerant carrier population. 

PACS numbers: 75.30.Ds, 75.10.Lp, 75.47.Lx 



I. INTRODUCTION 

Long-range ferromagnetic order mediated by inter- 
actions between itinerant and localized spins is es- 
tablished in many different materialsiii2i2^ The man- 
ganese oxides Ri_ x A x Mn03 (R=La,Pr,Nd,Sm,. . . and 
A = Ca, Ba, Sr, Pd, . . . ) are one prominent example.- 
The magnetic and transport properties of such itiner- 
ant magnetic systems are intimately related and, unlike 
in other ferromagnets, can potentially be controlled by 
tuning parameters like the itinerant carrier density. 

Learning how to control magnetization dynamics and 
relaxation is important for spintronics applications^^ 
One of the challenges facing future magnetic devices and 
memories concerns their speed, which is governed by the 
dynamics of the collective spin. For small deviations 
from equilibrium, the dynamical magnetic properties 
are determined by the spin susceptibilityj&i&Siiii Within 
the Random Phase approximation (RPA) £lL12£L12J£*i!i. 
which gives the spin susceptibility to 0(1/5*) (S is 
the local spin magnitude), magnetization relaxation 
arises from the interplay between the dephasing of 
the itinerant carrier spin and the magnetic exchange 
interaction i 8 ' 9 : 10 i 17 Spin relaxation also arises from in- 
elastic scattering processes, such as magnon scattering 
with charge excitations. 14 ' 18 i 19 Experimental probes of 
such effects include neutron scattering, ferromagnetic 
resonance, and ultrafast magneto-optical pump-probe 
spcctroscopyJiSai^-^^i2i^^i21 The interpretation of 
such experiments requires the development of many- 
body theories of spin dynamics and relaxation. 

Our goal in this paper is to develop a theory that de- 
scribes the local spin Green's function 

« »= -W) < [£tq(<). S Q (0)] >, (1) 

which determines the transverse spin susceptibility. In 
the above we introduced the collective spin operators 



S™ = l/VA 7 Ej S " e " iq Rj I «= x ,y, z , wh ere Sj describe 
spins localized on N lattice sites at positions Rj. S — 
S x HS y are the spin raising/lowering operators and (• • • ) 
denotes the grand canonical ensemble average. We con- 
sider a model Hamiltonian that accounts for the most 
important features common in a wide range of different 
itinerant ferromagnets, H = K + i/ C xch + -Haf + Hjj. 
K = X)kcr £ kcj co .Ck CT describes a band of itinerant car- 
riers, which in the manganites arises from the Mn d- 
states with e g symmetry, creates an electron with 
momentum k, spin a, and energy £k- To simplify the 
calculation of correlation effects common in many differ- 
ent physical systems, we consider a one-band model of 
n = 1 — x itinerant electrons per Mn atom that hop be- 
tween nearest-neighbour ltice sites. The electron concen- 
tration is described by the filling factor n = N e /N, where 
N e is the total number of electrons, which varies from to 
1. £k = — tjk, where 7k = 2J2i=i cos(kia). d is the sys- 
tem dimensionality and a the lattice constant (a = H = 1 
from now on). In the manganites, 0.5< n <0.8 in the 
metallic ferromagnetic regime of interest here. Our cal- 
culation can be extended to include the bandstructure of 
individual materials. 

In momentum space, the magnetic exchange interac- 
tion between the local and itinerant spins is given by 

k,qcr 

-^E{ s A-^i+^) ( 2 ) 

k,q 

where er = ±1. In the manganites, Eq.Q describes the 
Hund's rule onsite interaction between the e g carrier spin 
and the S = 3/2 local magnetic moment of the three 
electrons in the tightly bound ti g Mn orbitals. J ~ 2eV 
and 0.2eV < t < 0.5eV are typical values quoted in the 
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literature, which give 4 < J/t < 10. 29 The Hamiltonian 
K +i?exch defines the simple double-exchange modeh 28 i 29 
Here we add to the above minimal double exchange 
Hamiltonian two ubiquitous interactions; 

HjJ = N H C k+qT C L-ql C k'iC kT , (3) 
kk'q 

is the on-site (Hubbard) Coulomb repulsion. The 
Coulomb energy U ~ 3.5 — 8eV is the largest energy scale 
in the manganites^ 9 - This Hubbard interaction is gener- 
ally hard to treat and its effects on the spin dynamics 
have received less attention! 14 i 15 i 30 

Haf = Jaf y^Tk'S'k'S'-k 

k 

+ ^E7k(^ k + ^ k + ^ k ^ k ) (4) 

k 

is the direct super-exchange anti-ferromagnetic in- 
teraction between the nearest-neighbor local spins. 
Jaf ~0.01t is weak in the manganitesi^ 9 . 

Given the large values of J/ 1 in most systems of inter- 
est, many theories start from the strong coupling limit 
(J — ► oo) of the above Hamiltonian. In this limit, 
the itinerant carriers can hop on a site only if their 
spin is parallel to the local spin there. The kinetic en- 
ergy is reduced when all spins are parallel (double ex- 
change mechanism- 28 ), which favors the fully polarized 
half-metallic state \F). In the classical limit, S — > oo, the 
problem can be mapped to an effective nearest neighbor 
Heisenberg model with ferromagnetic interactions. The 
lowest order (0(1/ S)) quantum corrections are described 
by the RPA, which for strong couplings gives a disper- 
sion that again coincides with that of the nearest neigh- 
bor Heisenberg ferromagnet^ We may therefore assess 
the importance of correlations and quantum fluctuations 
beyond 0(1/ S) by fitting the Heisenberg dispersion to 
the experimental result and looking for deviations. 

For electron concentrations n > 0.7, initial measure- 
ments found nearest neighbor Heisenberg model spin 
dynamics. 31 However, later experiments reported devia- 
tions that increase strongly for n <0.7,2M2>Mi2£i2&i2Z>2Ma 
The Heisenberg model with nearest neigbor interactions, 
Ji , misses a pronounced softening near the Brillouin zone 
boundary. This softening is accompanied by a strong in- 
crease in the spin-wave damping as we approach the zone 
boundary. The experimental dispersion could be fitted 
by adding a fourth-nearest-neighbour ferromagnetic ex- 
change interaction, J4, to J\ while keeping J$= J-2=0. 39 

The above experimental observations reveal a new spin 
dynamics and nonlocal correlations that are not cap- 
tured by the strong coupling limit of the double ex- 
change model (which favors local correlations). Sev- 
eral scenarios have been put forward. The proposed 
mechanisms involve, among others, magnon scattering 
with orbital degrees of freedom ; 36 ! 40 Fermi sea pairs ; 30 ! 41 
or phonon o 33 ' 36 ! 42 , disorder effects^ 3 - bandstructure 



effects, 44 Hubbard interactions 1 ^ and correlations 3 - -, and 
the energetic overlap between spin-wave modes and the 
Stoner continuum^ The observed pronounced depen- 
dence of the spin-wave dynamics on the carrier con- 
centration puts stringent conditions on the theory. Ye 
et. al£^ argued that none of the mechanisms proposed so 
far can fully account for all aspects of this spin dynamics. 

The purpose of this paper is two-fold. First, we study 
the momentum dependence of the spin-wave dephasing 
rate, with the focus on the role of correlations between 
the spin and charge excitations and on their interplay 
with the itinerant spin dephasing. We demonstrate that 
both the magnitude and the momentum dependence of 
the spin relaxation rates depend sensitively on the carrier 
concentration and on correlations due to both J and U. 
We show that the on-site Coulomb repulsion U changes 
the momentum dependence of the spin-wave dephasing 
rate in a qualitative way that resembles the experimental 
results. We also show that our results depend on the in- 
terplay of U with the itinerant carrier spin lifetime, which 
is finite in some systems due to interactions not included 
in our Hamiltonian* 7 ^ We compare with the 1/S expan- 
sion and other approximations and find that three-body 
correlations between spin and electron-hole pair excita- 
tions play an important role. Finally, we show that the 
magnetization relaxation can be controlled by tuning the 
carrier density. We obtain changes in the spin relaxation 
with n that correlate with corresponding changes in the 
spin-wave softening and non-Heisenberg behavior. 

Second, we develop and test a general method for de- 
scribing spin dynamics. For this we use a truncation 
scheme of the infinite hierarchy of Green's function equa- 
tions of motion based on an expansion in terms of cor- 
relations. Our scheme treats the full dynamics induced 
by the correlations between a given number of elemen- 
tary excitations. Here we describe all correlations be- 
tween a local or carrier spin excitation and an electron- 
hole Fermi sea pair and obtain the exact solution within 
the sub-space of states with up to one Fermi sea excita- 
tion. Similar to Refs. |30l l 4lT |. our method becomes exact 
in the limits of one electron (N e =l, n=l/N), half fill- 
ing (N e —N, n—1), and in the atomic limit (t=0 for any 
n). It interpolates between the weak and strong coupling 
limits and agrees with exact diagonalization results for 
the spin-wave dispersion? 4 ^ Finally, it retains its varia- 
tional nature in the limit of zero relaxation rates, which 
provides a rigorous bound for the spin-wave softening 
and ferromagnetic phase boundary. Our approach, used 
before in the context of the Hubbard Hamiltonian ) 47 i 48 
is in the same spirit as the projection and factorization 
scheme of Ref.[49|, used to calculate the ultrafast non- 
linear optical response of systems with a strongly cor- 
related ground state, and the correlation expansion of 
Ref.[5(|. It may be extended to study spin correlations 
in non-equilibrium system o 17 ! 27 and the magnetization 
dynamics of (III,Mn)V semiconductors ^ 

The rest of this paper is organized as follows. In Sec- 
tion |H] we discuss the Green's function truncation scheme 
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and derive a closed system of equations of motion that 
determine the spin Green's function. In Section IIIII we 
obtain the spin self-energy and separate the RPA contri- 
bution from the contributions of correlations due to J and 
U. In Section HVl we discuss our numerical results for the 
spin-wave dephasing rate and dispersion and compare 
different approximations. In Section flV Al we consider the 
minimal double exchange model with U=Jaf=0, while 
in Section IIV Bl we study how U and Jaf change the 
picture. We end with our conclusions in Section fVl 



II. TRUNCATION OF GREEN'S FUNCTION 
HIERARCHY 

In this section we obtain the equations of motion that 
determine the Green's function Eq.(fT|) and the spin sus- 
ceptibility. The many-body interactions H exc h and Hu 
introduce an infinite hierarchy of coupled equations of 
motion that involve higher Green's functions of the form 

< A >= -i6(t) < [j4(£),Sq(0)] >, (5) 

where A(t) = exp (ziit)Aexp (—iHt) are many-body 
Heisenberg operators. To truncate this hierarchy, we ap- 
proximate the higher Green's functions by systematically 
adding correlations among any given number of elemen- 
tary excitations. To lowest order, the RPA describes un- 
corrected quasi-particles. At the next level, we include 
all correlations between any two elementary excitations, 
which determine the inelastic dephasing rate. 

Ref.[47] used a three-body scattering theory to calcu- 
late the electron Green's function of the Hubbard Hamil- 
tonian. In one dimension, the results obtained this way 
were in excellent agreement with the exact Bethe ansatz 
solution ! 47 ! 48 In Ref.[49[, a similar method was used to 
calculate the density matrix that describes the coher- 
ent ultrafast nonlinear optical dynamics of the Quantum 
Hall system. Refs. 52.53.54] calculated the Fermi Edge 
Singularity in doped semiconductor quantum wells us- 
ing an analogous approach. In this paper, we establish 
the correspondence with a factorization scheme of higher 
Green's functions. For simplicity we restrict to zero tem- 
perature, where Eq.([5]) involves the ground state average 
value. 

In the case of ferromagnetic exchange interaction as in 
the manganites, the fully polarized state 

i^Hn c >>(S)i^--->' (6) 

V 

is an exact eigenstate of the many-body Hamiltonian H . 
In Eq.©, |0) is the vacuum state and \S, S, ■ ■ ■ ) describes 
local spins with S z = S on all lattice sites. From now on, 
the indices fi, v, ■ ■ ■ denote states occupied in \F), while 
a, (3, ■ ■ ■ denote empty states. In the parameter range 
of interest here, \F) is the ground state i 30 ' 41 Using the 
properties H\F) = (we choose the eigenvalue of \F) as 
the zero of energy) and (-F^q = 0, both of which stem 



from the fact that \F) is the state with maximum spin, 
we obtain from Eq.© 

<$:A-»=~i8(t)<F\Ae~ lHt S Q \F> . (7) 

The Green's function <§; A ^> is then given by the ampli- 
tude of the time-evolved state SqI-F >■ The hierarchy 
of Green's function equations of motion is equivalent to 
solving the time-dependent Schrodinger equation. How- 
ever, Green's functions also treat dephasing and relax- 
ation, and can describe phenomenologically the effects 
of coupling to degrees of freedom not included in the 
Hamiltonian H by introducing phenomenological damp- 
ing rates. The coupling of -C A ^> to higher Green's 
functions is determined by the states HA\F). Truncation 
of the equations of motion hierarchy can be achieved by 
expanding HA\F) in a truncated basis, which gives the 
exact solution within a subspace of states. 

We start with the equation of motion for the spin 
Green's function Eq.(TT]), obtained after straightforward 
algebra by using Eq.JT]) and the properties of \F): 

(id t - ^ - u^j « » 

= 2SS(t) - ^= « c U c -+Ql » 

+4^ E « s i-»-Q c U c «i »< ( g ) 

OLV 

and 

c^ F = 2 Jaf5( 7 q - 7b) (9) 

is the spin-wave energy due to H af ■ The same result can 
alternatively be obtained by decomposing the Green's 
functions contributing to <C [S^q,H] ^> into correlated 
and uncorrelated parts after using the identity 

= < C L Q - qff C l-') « » 

+ (S$) « c{_ Q _^ a , » 

+ « AS£A[4_ Q _ q(TCW ] », (10) 

where S n are the components of the local spin and 
AA = A — (A) describes the quantum fluctuations of 
A. The Green's function < AS^A^cv] > is the cor- 
related part of S n c i a c a i 

The Green's function <§; cJ,*<v_|_q| 3> describes the 
itinerant carrier spin dynamics and satisfies the follow- 
ing equation of motion, obtained after straightforward 
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algebra by using Eq.([7|) and the properties of \F): 
(id t - £i,+q + e v - JS - nil) < cJ^c+qj, >= 

J ot u t 
-^7f « 5 -Q » - ]v E « <T c m+Q1 » 

V M 
J 



and 



2ViV 

-E 



E 



(11) 



where we defined the correlated part of the four-particle 
Green's function as 

•C efc2djC4 » C =<C cfc2c|c4 » 
~(cl c 2> < 4 C 4 > -(4 C 4) < c l c 2 > 

-(c 2 4) < c|c 4 > +<c|c 4 ) < 4 C 2 > • (I 2 ) 



Alternatively, Eq. (TlTj) can be derived by using Eq.JTI 
to decompose the Green's function <C [c' u ^c u+ qi, H] S>. 
The first line on the right hand side (rhs) of Eqs.© 
and (jTTJ) gives the RPA result, which neglects all non- 
factorizable (correlated) contributions to Eqs. lfTX))) and 
(fT2]) (Tyablikov approximation 5 - 5 -) . 

Eqs.© and (JTTJ) reduce the calculation of the spin 
Green's function to that of two higher Green's func- 
tions, <C A5 + A[ctc|] 3> and -C c|c|c|c| ^> c . In a sys- 
tem with a general ground state, two additional Green's 
functions, < A5 z A[c|c i ] > and < AS + A[c[ci] >, 
also couple and describe ground state correlations. 51 
However, these vanish here since, for the groud state 
Eq.©, AS z \F)=0 and c{cj.|F)=0. For the same rea- 
son, <C c j J T c ^ < M Cct T ^—-^ 4t Cj -4| c "T and <C 
s+c U c a1 >=< A5+A[c|; t c qT ] >. Also, from Eq.© 
we see that < A »= for any A such that (F |A = 0. 

The Green's function -C S' + cLc Q ,f S> describes the cor- 
relations between a magnon, an electron, and a Fermi sea 
hole (three-body correlations), while the Green's func- 
tion -c c^CQ+n+u-aic^Ca^ ^> c describes the correla- 
tions between two Fermi sea holes and two electrons 
of opposite spin. These two Green's functions are ob- 
tained from the following equations of motion, derived 
from Eq.([7]) after using the properties of \F): 



idt - e a +e v — 



J 

2N 
J 



rt 



JS 



Q 



N 



'Q+v-aj ^ S a _ v _QC ]/ ^C a -\ ^> 



^E <<: 4-^q4t c /5t» 



2N 
JS 



<C ^a-ju— Q C mT Cqi T > 



(13) 



J 



r + t 



J7 

' A 



^2 ^ c i'T c Q+M'+^-"i c lT c "T >c 



(14) 



Eqs. (fT3")) and (I14|) describe vertex corrections to the 
carrier-spin interaction. The first line on the rhs of 
Eq. (|13j) gives the Born approximation. The second and 
third lines describe vertex corrections due to the multi- 
ple scattering of the localized spin with the Fermi sea 
pair electron (second line) and hole (third line). Ne- 
glecting the third line corresponds to assuming a non- 
interacting (static* 5 - 2 -) Fermi sea, equivalent to summing 
only the electron-magnon ladder diagrams (two-body 
ladder approximation). — By also including the hole mul- 
tiple scattering processes (third line), we treat exactly 
all correlations between local spin, electron, and hole, a 
three-body problem. The last line on the rhs of Eq. (fT3|) 
comes from correlations between two electrons and two 
holes, described by Eq.(fT4|). We note that U introduces 
new correlations among all four of the above particles, 
described by the last four lines on the rhs of Eq. (fT4|) . 

To obtain the above closed system of equations, we 
neglected the coupling to Green's functions of the form 
<C Ac v fC' ^Caicpi where A = S + or c^cj.. These 
neglected Green's functions describe multi-particle cor- 
relations between two Fermi sea pairs and a local spin or 
carrier spin-flip excitation, which contribute to higher or- 
der in 1/S. Alternatively, we can arrive at the same result 
by decomposing the Green's functions -C S+cjcjcjC-f 3> 

and -c cjc|cjcjc|C| 3> into uncorrelated and correlated 
parts, by separating out all possible factorizable contri- 
butions similar to Ref.[5(|, and neglecting the fully cor- 
related contributions that describe correlations among 
three excitations. This correlation expansion neglects the 
contribution of states with two or more Fermi sea pair ex- 
citation and corresponds to an exact calculation of the 
Green's functions within the sub-space of states with up 
to one Fermi sea pair. As discussed e.g. in Refs. [47ll48j 
and implied by Eq.([7]), in the limit 7,T^0 the exact cal- 
culation of the Green's function within a given subspace 
it equivalent to the variational calculation of the spin- 
wave energy using a variational wavefunction that is a 
linear combination of the states that span the subspace 
(obtained in Refs. [3pK0 ] for the problem at hand). One 
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could include multipair correlations, e.g., by extending 
the approach of Refs. [SalHl. l54j |. 



III. SPIN SELF ENERGY 

The spin self energy can be calculated by solving the 
equations of motion derived above by Fourier transfor- 
mation. Eq.© gives 



25 



'-' " u> — Wq F — S(o>, Q) ' 



(15) 



where Q) is the self energy. Defining for convenience 



Xv(w, Q) 



G au (u>, Q) 



< 5; q > w 



(16) 
(17) 



and substituting into Eqs. (JHJ) , (fll|). (fT3| . and ((l4|) we 
express the self energy in the form 

^,Q) = ^-^E^+^E G - ( 19 ) 



where A = J 5 is the magnetic energy. We calculate S 
non-perturbatively in the interactions and 1/5 by solving 
the following coupled equations for X, G, and 

(nU + A + e v+ q - e„ - w) A„ - E X m 



v \ a / a// 



(20) 



Jn 



w + e v - e a - ~— - WQ F a+ „ + 17 ) G av 



2N 



1 + E G /3" ~~ E ^ 



(21) 



where 7^0, and 



(nC/ + A + £ a + e Q+f<+ „_ a - e,, - e„ - w) $^ 



u 
'n 



(22) 



First we consider the RPA self energy, obtained by 
setting G = $ = in the above equations. We can then 
solve Eq. ([20|) analytically after noting that its solution 
has the form 



X 



RPA 



(w,Q) 



(23) 



nU + A + e u+ q — e v — uj 
Substituting the above expression into Eq. (j2"0|) we obtain 
J 1 



X 



RPA 



which gives after some straightforward algebra 
J 1 



(24) 



2y/N A + U„q + e u+ d -e v -ui 
where we introduced the Coulomb-induced energy 



, (25) 



N ^"nlI + A-iT + e M+Q 



(26) 



Substituting Eq.([2g]) into Eq.([l9"l) after setting G = 0we 
obtain the RPA self energy 



^RPA, 



w,Q) 



J 

2N 



E 



■ CO 



A-iT + U u 



Q 



-H-Q 



_ (27) 
Following Ref. [7[ , in the above equations we added a phe- 
nomenological relaxation rate T that describes the itin- 
erant spin lifetime due to interactions not included in 
our Hamiltonian H. This result can alternatively be ob- 
tained by substituting A by A — iT, as derived with the 
Lindblad semigroup method in Ref. [l7T] . In the intrinsic 
system described by the Hamiltonian H , T^O. 

We now turn to the self-energy due to the cor- 
relations. By formally solving Eq. ([20|) for X v and 
substituting into Eq. (fP9|) . we separate the self-energy 
into RPA and correlated contributions, £(o>, Q) = 
S RPA (w, Q) + E corr (w, Q). After some algebra we obtain 
that £ corr (w, Q) = S c 7 orr (w, Q) + S^ orr (w, Q). 



V^corr \ A s~i 



u vCl 



— e v — uj 



iT + U v q + £ih-q - e v 



■ CO 



-(28) 



is the contribution of the Fermi sea-magnon correlations 
due to J, described by the Green's function G, Eq. (|2"Tj) . 

[/ v^.„ A - iT 



TV 3 / 2 



A — iT + U V C, + £„+Q — Ev - W 



(29) 



is the contribution of the Fermi sea pair-carrier spin-flip 
four-particle correlations described by and discussed 
in Ref.[30(. This latter contribution vanishes for U=0. 



IV. NUMERICAL RESULTS 

In this section we present our numerical results for the 
spin-wave dispersion, which we obtain by solving self- 
consistently the equation 



af 



Re£( WQ ,Q), 



(30) 
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FIG. 1: (Color online) Spin-wave dispersion for different val- 
ues of y/t for J — St,n — 0.6,17 = Jaf — T — 0. Inset: 
Dispersion along the F — > X direction. 



and the spin-wave dephasing rate Tq, determined by 
ImE(wQ,Q). We focus on the inelastic contribution to 
the dephasing rate, due to the scattering between spin 
and charge excitations: Tq 



Tq, where 



1 Q 



-Im£ c j or >Q,Q) 



(31) 



is the contribution due to the magnon-Fermi sea pair 
correlations described by G, Eq.(|2"8j). and 



(32) 



r^ = -imsr r (^Q,Q) 



is the contribution due to the carrier spin flip-Fermi sea 
pair correlations described by <£>, Eq. (j2"2"l) . An additional 
elastic contribution to the spin-wave lifetime can come 
from the imaginary part of the RPA self energy, which 
however vanishes in the limit r — ► due to the finite car- 
rier spin-flip excitation energy. Within the RPA, spin- 
wave dephasing can only arise from the interplay between 
the magnetic exchange interaction and carrier spin de- 
phasing via external couplings 

Below we study the momentum dependence of Tq 
along different directions in the Brillouin zone. We focus, 
in particular, on the directions T-X, T-M, and X-M, 
where T = (0, 0), X = (it, 0), and M = (tt, it). Our calcu- 
lations were performed on a 20x20 square lattice, which 
as shown in Refs.[4l| and [3Cj gives good convergence 
to the thermodynamic limit. Any small size effects are 
washed out when the relaxation rate 7 in Eq. (f2~I| exceeds 
the energy spacing. 



A. Minimal double— exchange model 

First we consider the simple double exchange Hamil- 
tonian and set U = Jaf = T = 0. FigQ] shows the spin- 
wave dispersion loq as function of the damping rate 7 in 



Eq. (f2T]) . With decreasing 7, our results converge to the 
spin-wave dispersion obtained variationally in Ref.(4l| 
and the 7 — * limit. Similar to the experiment^ our 
calculation can then be fitted to the dispersion of the 
Heisenberg Hamiltonian with first- and fourth-nearest- 
neighbor spin interactions J\ and For intermedi- 
ate concentrations, our calculation gives a pronounced 
spin- wave softening at the X point, described by J4, as 
compared to both the RPA and to the fit to the nearest 
neighbor Heisenberg dispersion. The effects of a finite 7 
are most pronounced along the direction r — > X: with 
increasing 7, the time evolution described by the Green's 
function G, which determines the spin-wave softening, 
is suppressed and thus the dispersion starts to approach 
the RPA (G = 0) result. From now on we fix 7 = 0.2t, 
which as can be seen in FigQ] is close to the 7^0 limit. 

Fig[2] demonstrates the important role of correlations 
due to spin-charge interactions on both the spin-wave 
dispersion and dephasing rate. The spin-wave energies 
and lifetimes differ markedly depending on the approxi- 
mation used to treat the correlations. The latter deter- 
mine the differences from the RPA, which describes non- 
interacting spin-waves (G = $=0). The RPA grossly 
underestimates the softening and does not give any spin 
damping in the limit T — > 0. 

By neglecting the Green's function <!>, we obtain spin- 
wave energies closer to the RPA (see Fig|3Ja)). As seen in 
Figj^b), this approximation, which only treats the scat- 
tering of local spins with the Fermi sea, gives a very small 
damping rate. On the other hand, the non-variational 
0(1/ S 2 ) approximation discussed in Refs. |l4lll9l |. which 
treats magnon-Fermi sea scattering within the Born ap- 
proximation and Fermi's golden rule, strongly overesti- 
mates the softening, while at the same time predicting 
only a small damping rate (see Fig[2]). 

By adding to the 0(1/ S 2 ) result the effects of the mul- 
tiple scattering of the magnon with the Fermi sea pair 
electron, while still neglecting the magnon-Fermi sea hole 
interactions, we obtain a non-variational two-body ap- 
proximation of the vertex corrections equivalent to sum- 
ming the magnon-electron ladder diagrams. 47 As can be 
seen in FigEl this ladder approximation gives very large 
softening and damping, much larger than the predictions 
of the full calculation (which is variational in the limit 
7,T^0 considered here). The latter treats, in addition 
to the magnon-electron interactions, the multiple scat- 
tering of the Fermi sea pair hole with the magnon. The 
large differences between the ladder and full calculation 
results demonstrate the importance of three-body corre- 
lations between magnon, electron, and hole in the param- 
eter regime of interest in the manganites. We conclude 
based on Fig[2] that all correlations between spin, elec- 
tron, and hole must be treated on an equal basis. The 
variational nature of our full calculation of the spin-wave 
energies in the limit 7, T— >0 has the advantage of provid- 
ing a rigorous limit of the magntitude of the softening, 
unlike for the ladder or 1/S expansion results. 

Fig[3] shows the behavior of ojq and Tq for different 
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FIG. 2: (Color online) Comparison of the different approx- 
imations for treating the correlations for (7=0. (a) Spin- 
wave dispersion, (b) Inelastic spin-wave dephasing rate, n — 
0.6, J = 8t, 7 = 0.2i, U = Jaf = T = 0. 



FIG. 3: (Color online) Dependence of the spin dynamics 
on the interaction strength J. Spin-wave dispersion (a) and 
dephasing rate (b). n = 0.6, 7 = 0.2t, r = U = Jaf = 0. 



interactions J. With increasing interaction strength, the 
spin-wave energies increase and the ferromagnetic phase 
becomes more stable. This hardening with J is accom- 
panied by a corresponding increase in the spin-wave life- 
time. The above changes are stronger along the T-X 
direction, where the non-Heisenberg behavior and soft- 
ening are pronounced, and depend nonlinearly on J/t. 
The overall momentum dependence, however, remains 
the same for all J. 

To interpret the above results, we turn to the Green's 
function equations of motion and note that, for U=0, 
E u =0. After solving Eqs.pO]) and J32J for X and $ and 
substituting into Eq. (|2Tj) we obtain that 



^auG aL/ 



2N A - iT + e u+Q -s v -u 



EQ+H+u- 



A — iT + SQ+ft+u-a — £fi + e a — e u — 



(33) 



where 



VL au — lu + i-y - (e a - e v ) 



J 

~2N 



E 



£QH-/x+i/— a ^ f-i 1 £q 



A - iT + £ Q 



+fj,+v— a 



- £„ 



43.4) 



£ n — £ v —LU 



and 7,T ->0. 

Spin-wave dephasing results from the scattering of the 
magnon of momentum Q to momentum Q + v — a while 



an electron is excited from the state v inside the Fermi sea 
to the empty state a. In the limit 7^0, this magnon- 
Fermi sea scattering process must satisfy the energy con- 
servation condition f2 Ql/ = 0, i.e. the initial magnon en- 
ergy, to — luq, must equal the final state energy that in- 
cludes the Fermi sea pair energy and the Q+is— a magnon 
energy. The final state magnon energy, given by the last 
term in Eq. (f34|) . comes from the coupling of G to <I>. For 
$=0, this spin-wave energy is replaced by the local spin 
excitation energy Jn/2. The scattering of small energy 
Fermi sea pair excitations from right below to right above 
the Fermi surface dominates the spin-wave lifetime. The 
density of states and characteristic momenta of such pair 
excitations depends on the shape of the Fermi surface 
and therefore on the carrier concentration. 

To derive the 0(l/5 3 ) dephasing rate t 14 ' 19 we neglect 
all rescattering terms on the rhs of Eci. (f33|) (oc G). Sub- 
stituting the expression for G obtained this way into 
Eq. (|28p , we obtain the Born approximation self energy 



— 



J 2 
4N 2 



E 



— Ev — U) 



if- (35) 



Expanding in terms of 1/S while keeping A = JS 
fixed and substituting u = Wq , where Wq' denotes the 
0(1/5*) spin-wave energy, we obtain the lowest order con- 
tribution to the self-energy imaginary part: 
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FIG. 4: (Color online) Dependence of the spin-wave dephasing rate on the carrier concentration n for (7 = 0: comparison 
between the different approximations of the correlations J — 8t, 7 = 0.2i,F = U = Jaf = 0. 



where 

= ~ ( £ a -e u + u£>_ a+l /J - ij. (37) 

The above result corresponds to the Fermi's golden rule 
description of the magnon lifetime. Its large difference 
from our full calculation, demonstrated by Fig(2l is due 
to the magnon electron and magnon-hole multiple in- 
teractions (vertex corrections), described by the terms 
proportional to G on the rhs of Eq. (|33|) . The compari- 
son between the different approximations shows that, in 
the parameter regime of interest in the manganites, the 
vertex corrections due to three-body correlations renor- 
malize significantly the magnon-carrier scattering. 

We finally turn to the dependence of the spin relax- 
ation on the carrier concentration. Fig|4ja) demonstrates 
a strong n-dependence of Tq, which correlates with an 
analogous dependence of wq and the spin-wave soften- 
ing discussed in Ref.[3(|. As n decreases and the soft- 
ening (non-Heisenberg behavior) disappears, the spin- 
wave lifetime increases while its momentum dependence 
changes. For intermediate n, Tq displays two sharp 
peaks and a dip as function of momentum. For small 



n, the overall Tq decreases and the positions of its max- 
ima and minima change. 

To see this concentration dependence in more detail, 
we note that, for n — 0.72, the spin-wave damping is 
maximized for Q ~ (tt,tt/2), between X and M, and 
Q ~ (tt, 7r/2), between M and T, while it is minimized 
close to the M point. As n decreases to intermedi- 
ate values, the first of the above maxima approaches 
the X-point while the second maximum shifts closer to 
Q = (tt, 7t/2). For smaller n, the dip close to the M point 
turns into a maximum. As a result of this n-dependence, 
Tq/luq becomes quite large in the direction r — * X 
for intermediate n, which implies that these spin-wave 
quasi-particles interact strongly with the Fermi sea. For 
small n, Tq/wq decreases again. The changes in the mo- 
mentum dependence with n are related to the changes in 
the shape and position of the Fermi surface, which is lo- 
cated close to the Brillouin zone boundary for the higher 
n but moves towards the center of the Brillouin zone as 
n decreases. As a result, the phase space available for 
magnon-carrier scattering changes drastically with n. 

FigU also compares the concentration dependence of 
Tq predicted by the different approximations of the spin- 
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FIG. 5: (Color online) Role of the super-exchange interaction 
Jaf. (a) Spin-wave dispersion and (b) Spin-wave damping 
rate. J = 2t, U = 10*, n = 0.6, T = 0M, 7 = 0.2i 



charge interactions. By comparing the full three-body 
calculation with the 0(l/5 3 ) Fermi's golden rule result, 
it is clear that the spin-wave damping is grossly underes- 
timated by the perturbative 1/S expansion for all concen- 
trations. FigslUb), obtained by setting <I>=0, fails com- 
pletely to capture the correct concentration dependence 
(compare Figsllja) and0|b)). It also predicts very small 
dephasing rates for all n. The above approximation ne- 
glects the interactions between Fermi sea pair and carrier 
spin-flip excitations. We therefore conclude that such 
carrier-carrier interactions strongly affect the magneti- 
zation relaxation. Finally, the comparison of FigsHJa) 
and (c) shows that the two-body ladder approximation 
grossly overestimates the spin-wave damping for inter- 
mediate or high n, while the discrepancies from the full 
three-body calculation decrease for small n. We con- 
clude based on Figf4] that the collective spin relaxation 
predicted by the minimal double-exchange model can be 
controlled by tuning the carrier concentration n, by dop- 
ing or with external probes. Such tuning is heavily influ- 
enced by the correlations, which must be treated accu- 
rately in order to capture even the correct order of mag- 
nitude and momentum dependence of the spin dephasing 
rate for all concentrations. 



B. The role of the Coulomb Repulsion 

In this section we study how the on-site Coulomb re- 
pulsion, U, and direct superexchange interaction, Jaf, 
affect the spin-wave energies and lifetimes. FigJS] com- 



pares the results obtained for different values of Jaf/£ 
within the range < Jaf < 0.012t relevant to the 
manganites. 29 Jaf leads to an overall softening of the 
spin-wave energies. These eventually turn negative, im- 
plying instability of the fully polarized ferromagnetic 
phase. However, Jaf preserves the nearest neighbor 
Heisenberg model momentum dependence, unlike for the 
softening observed in the experiment^ Also, it does not 
affect the spin damping in a significant way. We take 
Jaf = from now on. 

As demonstrated by FigdJJ the effects of the on-site 
Coulomb (Hubbard) repulsion U are more significant. 
FigEIa) shows the dependence of the spin-wave disper- 
sion on U. With increasing U, the spin-wave softening 
and deviations from the Heisenberg model dispersion di- 
minish as the ferromagnetic phase becomes more stable. 
As discussed in Ref.[3(|, our calculated dispersion can 
be fitted to the Heisenberg model dispersion with first- 
and fourth-nearest-neighbor interactions, similar to the 
experiment. 39 FigJHJb) demonstrates qualitative changes 
in the overall momentum dependence of the spin-wave 
damping as compared to the minimal double exchange 
model. In particular, if the carrier spin is conserved 
(r=0), Tq is maximum at the M point, while the damp- 
ing at the X point is smaller. In contrast, for U = and 
intermediate concentrations, the dephasing rate displays 
a dip at the M point and is maximum close to the X 
point (see e.g. FigJSJ. The double-peak momentum de- 
pendence of Tq for U=0 can be recovered for large U only 
by introducing a sufficiently large itinerant spin damping 
rate T (see FigJ^c)). The experiment of Ref . [39j] observed 
an increase in the spin-wave damping along T — > M that 
exceeds the corresponding increase along r — > X , similar 
to our results for finite U and T=0 (Figjf^b)). Our calcu- 
lations show that such behavior of the spin damping can 
be attributed to the Coulomb repulsion U in the intrinsic 
system (r,7 — >0) described by the Hamiltonian H. 

FigO also demonstrates a qualitative difference in the 
dependence of Tq on U between the cases of small 
and large itinerant spin damping. For r=0, the spin- 
wave dephasing rate increases with U (FigjBJb)), while 
for r=0.5£ it decreases with U close to the A-point 
(FigJBfc)). This result indicates that the magnetiza- 
tion relaxation may depend on the interplay between 
Coulomb repulsion and the dephasing of the itinerant 
spin via the coupling to an external bath. 

To interpret the above behaviors, we plot in Fig[7]the 
two contributions to Tq obtained from the self energies 
Eqs. pS)) and (|2"9"j) for zero and finite T. T J is determined 
by G, while T u is determined by For r=0, FigH^a) 
shows that T J and T u are comparable in magnitude, 
since in this case they both arise from ImG. On the other 
hand, as T increases, the relative magnitude of T J and 
T u changes and the latter dominates (see FigfTfb)). This 
enhancement of T u arises from the additional contribu- 
tion to Im<i>, Ea.(|22p. obtained by adding the relaxation 
rate T to Eq. ([!?!?)) . Even though T J continues to have the 
same momentum dependence as for r=0, T u does not 
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FIG. 6: (Color online) Effects of the on-site Coulomb re- 
pulsion U. (a) Spin-wave dispersion, (b) Spin-wave damp- 
ing rate for F = 0, and (c) (b) Spin-wave damping rate for 
r=0.5t. n = 0.6, J = 2t, 7 = 0.2t, Jaf = 0. 



(compare FigsE^a) and[7jb)). 

Fig[8] shows the dependence of the two contributions 
to the inelastic dephasing rate on the Coulomb repulsion 
for large U and different T. T J increases and then satu- 
rates with increasing U. For T=0 (intrinsic system), T u 
increases with U and eventually exceeds T J . Unlike for 
T J , the dependence of T u on the momentum and on U is 
qualitatively different for large and small T. For exam- 
ple, T u decreases with U at the X-point for large F but 
increases for T=0. The above behavior of T u dominates 
the total dephasing rate for large T. 

FiglH a ) shows the transition in the momentum de- 
pendence of Tq as U increases for r=0. This transition 
occurs around U ~6t, where the double peak momen- 
tum dependence for [7=0, with a dip at the M-point, 
changes into a peak at the M-point. As can be seen in 
FigJSJb), the above transition arises from the changes in 
the behavior of T J , determined by the Green's function 
G (Eq.flSTJ), that are introduced by the Coulomb repul- 
sion. Fig[5Jc), on the other hand, shows that the overall 
momentum dependence of T u remains approximately the 
same for all U. 

Finally, we turn to the possibility of controlling the 
magnetization relaxation by tuning the carrier concentra- 
tion n and study how the Hubbard repulsion U changes 
the picture as compared to the prediction of Fig0J FigfTUl 
shows the dependence of the dispersion and spin damp- 
ing rate on n within a wide range of concentrations, 
n = 0.7 — 0.1. As can be seen in FigfTOTa). the pro- 
nounced softening along the T—X direction disappears 
rapidly with decreasing n (or increasing hole doping x=l- 
n). For small values of n, the overall energies decrease 
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FIG. 7: (Color online) Contributions Fq and Tq to the in- 
elastic spin-wave dephasing rate Fq. (a) r = 0, (b) F=0.5t. 
n=0.6, J=2t, U=10t, j=0.2t, Jaf = 0. 



and the overall shape of the dispersion changes. 

The above n-dependence of the spin-wave dispersion 
correlates with corresponding changes in the dephas- 
ing rate. FigfTUTb). obtained for T=0, shows that the 
Coulomb repulsion U changes drastically the dependence 
of Tq on n as compared to the predictions of the mini- 
mal double exchange model, FigH^a). In this case, the 
spin dephasing rate increases as n decreases to intermedi- 
ate values, in a way correlated with the disappearance of 
the spin-wave softening and non-Heisenberg behavior. 
For small n, where the softening has disappeared and 
the overall energies start to decrease, the spin damping 
rate also decreases. Furthermore, unlike for U=0 or for 
large T, Tq displays a strong maximum at point M for 
all concentrations and is always weaker along the T—X 
direction. As can be seen by comparing FigsfTUFa) and 
(b), a sufficiently short itinerant spin lifetime (large T) 
changes drastically the momentum dependence of Tq and 
its dependence on n. 

We conclude based on FigdHlthat the magnitude and 
momentum dependence of Tq, as well as the spin-wave 
softening, can be controlled by tuning the carrier concen- 
tration n, via hole doping or by external means such as 
photoexcitation,iL2aiSi^^i2i^^i22 It is clear that the 
on-site Coulomb (Hubbard) repulsion plays a dominant 
role, by inducing new correlations and dynamics absent 
in the simple double-exchange model. Such correlations 
change drastically the momentum dependence and mag- 
nitude of Tq with varying n and must be treated in a 
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FIG. 8: (Color online) Dependence of Fq and Fq on Coulomb 
repulsion for large U. n=0.6, J=2t, 7=0.2t, Jaf = 0. 




FIG. 9: (Color online) Dependence of (a) Tq, (b) T q , (c) 
on the Coulomb repulsion for small U. n—0.6, J—4t, 7=0. 2t, 
Jaf = r = 0. 



consistent way in order to arrive at trustworthy conclu- 
sions and comparisons to experiment. Our results sug- 
gest that, as a first step, a systematic experimental study 
of the magnetization dynamics as function of doping x 
and a comparison to the theory is necessary in order to 
decide which many-body mechanisms dominate the col- 
lective magnetization dynamics and relaxation and learn 
how to control this dynamics for potential magnetic de- 
vice and spintronics applications. 



V. CONCLUSIONS 

In this paper we presented a general method for de- 
scribing the spin-wave dynamics and relaxation in itin- 
erant ferromagnets. This method is based on a correla- 
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FIG. 10: (Color online) Dependence of spin-wave dispersion 
and dephasing rate on the carrier concentration n for U=10t. 
(a) dispersion, (b) dephasing rate for F=0, (c) dephasing rate 
for r=0.5i. J = 2t, Jaf = 0, 7 = 0.2t 



tion expansion of the Green's function equations of mo- 
tion that systematically treats all correlations between 
any given number of elementary excitations. Using this 
method, we derived a closed system of equations that 
treat the magnetic exchange and Coulomb interactions 
non-perturbatively and solved it to obtain the Green's 
function that determines the transverse spin susceptibil- 
ity. Our results for the spin-wave dispersion reproduce 
previous variationa l 30 ^ 1 and exact diagonalizatio n 46 i 47 
results (in the limit T,j — >0) and therefore allow us to 
draw definite conclusions regarding the magnitude of the 
spin-wave softening. Using the properties of the fully po- 
larized Hartree-Fock ground state with maximum spin, 
we showed that our method gives the exact spin Green's 
function within a subspace of states that include up to 
one Fermi sea pair excitation. Our factorization scheme 
of the higher Green's functions also applies to other 
ground states. Our results recover the 1/S expansion re- 
sults as special case. We showed that, in the parameter 
regime of interest in the manganites, the latter approx- 
imation overestimates the spin-wave softening, while at 
the same time it grossly underestimates the spin-wave 
damping rate. Furthermore, by comparing with the lad- 
der approximation treatment of the vertex corrections to 
the magnon-carrier scattering, which treats the multiple 
magnon-electron scatterings while neglecting the interac- 
tions with a Fermi sea hole, we showed that three-body 
correlations between the magnon and an electron-hole 
Fermi sea pair excitation have an important effect on the 
spin relaxation. 

Using the above many-body theory, we calculated the 
inelastic spin-wave dephasing rate non-perturbatively 
in the interactions and 1/S (i.e. beyond the standard 
Fermi's Golden Rule). We showed that correlations be- 
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tween a carrier spin-flip excitation and a Fermi sea pair 
induced by the Coulomb repulsion U play a very impor- 
tant role in the parameter regime relevant to the man- 
ganites. We also showed that both the magnitude and 
momentum dependence of the spin-wave dephasing rate 
depend sensitively on the itinerant carrier concentration. 
This result implies the possibility of controlling the mag- 
netization relaxation in itinerant ferromagnets by tuning 
the carrier concentration, either via doping or by exter- 
nal means (e.g. photoexcitation or by using electric fields 
and currents or gates). We also argued that the interplay 
between on-site Coulomb (Hubbard) interaction and a fi- 
nite itinerant carrier spin lifetime, due to the coupling of 
the different spin states induced by spin-orbit or other 
interactions not included in our Hamiltonian, can affect 
our results in the realistic system. The momentum de- 
pendence of the spin-wave dephasing rate observed in 
recent experiments^ is consistent with the results of our 



calculation only for sufficiently large U and V — >0. In all 
other cases, we obtain a distinctly different double peak 
momentum dependence. Although the bandstructure of 
the relevant materials must be included in order to ar- 
rive at quantitative comparisons with the experiment, 
our calculation already demonstrates the crucial role of 
correlations. The agreement of the main trends in spin- 
wave softening and damping rate changes as function of n 
between our theory and the experiment suggests that the 
simple one-band model already contains the main inelas- 
tic scattering processes and correlations. Our results sug- 
gest that ultrafast magneto-optical pump-probe spec- 
troscopy experiments, which directly probe the changes 
in spin relaxation and dynamics induced by photoexcited 
carriers, 17 may provide new insight into the physics of the 
manganites and other itinerant ferromagnetic systems. 

This work was supported by the EU STREP program 
HYSWITCH. 
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